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Processing Electrophysiological Signals 
for the Monitoring of Alertness 

Abstract 

This research project is concerned with the use of mathematical 
techniques for processing EEG signals associated with varying states of 
alertness. It requires the development and implementation of advanced 
signal modeling and data processing techniques; especially designed 
for the representation and prediction of bioelectric signals useful 
as estimators of states of alertness. In particular, our goal is to 
develop and implement efficient techniques for processing and modeling 
of EEG signals to extract the characteristics of signals associated with 
varying states of alertness. New represent at -i nn« fov signals 

which will enhance the features in estimating the states of alertness 
are sought after^ Fast algorithjns for implementing real-time computations 
of alertness estimates have also been developed. In this report, we 
present a new realization of the phase-distortionless digital filter 
which approaches real-time filtering and a new transform for EEG signals. 
This transform not only provides new information for the alertness 
estimates but also can be performed in real time. We are also developing 
a statistical test for stationarity in EEG signals. This test provides 
a method for determining the duration of the EEG signals necessary in 
estimating the short-time power or energy spectra for nonstationary 
analysis of EEG signals. It also helps in extracting the dynamic 
properties of the signal process. 


iii 


I 



Processing: Electrophysiological Signals 
for the Monitoring of Alertness 

I . Introduction 

The goal of this project is to develop and implement advanced 
techniques for processing and modeling of EEG signals to extract the 
characteristics of signals associated with varying states of alertness. 
These techniques are especially designed for the representation and pre- 
diction of electrophj^siological signals useful as estimators of states 
of alertness. The purpose of signal representation is to transform 
the data into a domain which is more convenient, for some reason, to 
the particular user or users or for further analysis; e.g., data reduction 
or the enhancement of alertness estimates may be one of these goals. 

Fourier series or transform has proven useful in the analysis of 
EEG signals in that it enhances the characteristics of the rhythmatic 
activities in EEG signals. With the advent of f ast-Fourier-transform 
algorithm, the use of the Fourier spectral analysis has been widespread. 
However, it requires that the signal be at least wide-sense stationary. 
Since the EEG signals manifests the subject’s brain states, the statistical 
characteristics vary as the states of the subject change during the 
course of time. It is more natural to treat EEG signals as nonstationary 
processes as proposed or as treated in [l] [2], The general approach is 
to use a short-time power or energy spectral analysis. The success of 
this approach depends on the signal process remaining stationary during 
the duration for which the spectrum is obtained. The quantitative 
determination of this duration is important for both the accuracy and 
the efficiency of this analysis. We have used a sensitive test which 
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will not only determine the time interval of the signal in which 

significant change in certain statistical characteristics is unlikely 

I 

to occur, but also detect any rapid changes in the interval if these 
statistical characteristics do vary. An algorithm to implement this 
test on a computer has been developed. In addition, this test will 
also detect the transition of the brain states and aid in extracting 
the dynamic properties of the process. We shall describe this test 
and its implementation and some preliminary results obtained by applying 
it to the EEG data furnished by NASA-Ames Research Center, 

As proposed, we have also searched for other representations 
(transforms) of EEG signals which may furnish advantages over the Fourier 
transform. The Fourier spectrum is a signal representation in which 
the sinusoidal waveforms are used as basis functions. Different sets 
of basis functions will result in different spectra. The search for a 
new signal representation (transform) is tantamount to the quest for a 
new set of basis functions. In particular, we have used a set of 
orthonormalized square waves as basis functions to represent EEG signals. 
This set of square wavefoi'ms is known as the Walsh functions [3] . The 
resultant representation by using this set of square waveforms is thus 
called the Walsh transform. The Walsh transform is especially suitable 
for digital or digitized signals. The algorithm for Walsh transform is 
far faster than the f ast-Fourier-transform, and it can be easily imple- 
mented by hardwire programming. Real-time computation of the Walsh 
spectra is a cinch. From our analysis of the EEG data furnished by 
NASA-Ames Research Center, the Walsh spectra enhance certain features 
of the EEG signals. We shall describe our analysis of EEG signals by 
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the use of Walsh transform and demonstrate its advantages. 

In the phase-uncertainty model of EEG signals [4], there is the 
characteristic parameter b which is the propagation speed of the 
variance used as a measure of the uncertainty of the phase. We believe 
that this parameter b will be a good measure for alertness. In order 
to obtain a time-varying b from the EEG signals, we have used a 
sliding window and an estimator method rather than the correlation 
technique. The ultimate goal is to find a scalar relationship between 
the time-varying h and behavioral alertness measurements. It is 
essential to have a linear-phase filter so that the phase information 
in EEG will not be distorted after filtering. The previously used time- 
domain transversal filter does serve the purpose; however, it is 
extremely slow. We have implemented a frequency-domain linear-phase 
filter which operates in real time. The details of this filter will be 
described. 
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II. A Test for Stationarity in EEG Signals 

To treat EEG signals as nonstationary processes, many researchers 
have used the short-time spectral analysis. This technique depends on 
the validity of the assumption that the signal remains at least wide^ 
sense stationary in the interval when the short-time power density 
spectrum is being computed. For a process to be wide-sense stationary, 
its mean must be a constant and its autocorrelation function must be 
a function of the lagging interval alone. We have devised a sensitive 
method for detecting any changes of the stationarity conditions in the 
signal process for any intervals of time. The approach is to subject 
the mean and mean square of the signal process to randomness statistical 
tests in order to identify non-random variations of the parameters of 
interest in the course of time. This method will not only determine the 
epoch within which the signal process remains stationary but also will 
detect the transition of the brain states and thus will aid in extracting 
the dynamic properties of the process. 

The mean, variance, and mean square are calculated for each successive 
one-second segment of data. These parameters then form a set of arrays 
consisting of from 5 to 25 elements. Each array then is composed of 
statistics representing an interval of data from 5 to 25 seconds, all 
intervals starting at the same reference as shown below: 

1 2 3 4 5 

M . 1 I 1 1 

1 2 3 4 5 6 

.1 1 1 „l 1 1 1 

1 2 3 4 5 6 7 


1 2 3 4 5 6 
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Each of these arrays is then subjected to four tests of nonrandomness, 

since each test is sensitive to a different type of trend. The four tests 

! 

are : ■ 

a) Run test : 

b) Trend test 

c) Mean Square Successive Differences test, and 

d) Serial Correlation test. 

Each test assigns a parameter value to each array. The parameter values 
are then plotted versus interval length, and trends are displayed. If 
there is a nonstationary trend, it will manifest itself by a significant 
change in one or more of the test parameters. The procedure is block 
diagrammed in Figure 1. 
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Figure 1 Block diagram for tests of nonstationarity 

As an illustration, we have used the EEG data recorded while the 
human subjects with eyes closed were stimulated by stroboscopic flashes 
at the rate of 10 flashes per second for 50 seconds, then no stimulation 
for 50 seconds, then another 50 seconds of stimulation, etc. The EEG 
data were furnished by Dr. J. Anliker of NASA-Ames Research Center. 

In this case tlie Mean Square Successive Differences test was the 
most sensitive to the state of strobe--on or off. It showed a change 
in test parameter starting shortly after the change in strobe status, 
and lasting several seconds. This corresponds to expected latency 
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periods for photically stimulated subjects. 


Results 

Figures 2-8 illustrate the effect of the strobe turning off. At 
the point marked T on the horizontal axis (e.g., 25 seconds in Figure 
2), the strobe went off. This means (for Figure 2) the strobe was on 
for 25 Seconds following the reference time. In Figure 3 the reference 
time in the EEG data is 6.3 seconds later than in Figure 2, and, therefore, 


the T indicating strobe off occurs after 18.7 seconds (25 - 6.3). 

In other words, the different figures represent the results of sliding 
a 25-second window along the EEG data with 6.3-second increments. 

The dependent variable is the result of the Mean Square Successive 

2 

Differences test run on the mean (m), variance ) and mean square 
2 

( 4 ' ) for each interval from 5 to 25 seconds. Thus, in Figure 3 the first 
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last seven are combinations of pre- and post-strobed terminations . 

As can be seen, there is no trend before T in all figures (which 
is reasonable since all the points represent strobed data). Starting 
with Figure 3 a trend can be seen beginning approximately one second 
after the cessation of strobing, and continuing for approximately 6 
or 7 seconds. This is best seen in Figure 6. Fig^ure 5 is a repeat of 
Figure 4 with polynomial regression lines (4th order) fit to the data. 
Examining this series of figures we. observe that the effect of turning 
off the strobe is detected, and the nonstationary interval following this 
event is identified. 

Figures 7 and 8 show the results of MvSSD test for stimulus -off EEG 
signals where no trends have been detected. 
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Figure 6 MSSD test for EEG signals in 
the stimulus-on state and then the 
stimulus-off state 







III. 


Representations of EEG Signals 


In our search for new representations of EEG signals, we aimed 
at the efficiency of the representations and their enhancement of various 
desired features of the signal. The Karhunen-Loeve expansion offers 
the most efficient representation; however, it suffers from its difficulty 
in implementation, Fourier transform has proven to be useful in some 
EEG analysis. In the course of our work, we have developed and imple- 
mented a statistical representation in terms of amplitude, frequency, 
and phase of EEG signals [s] and the representation in terms of in- 
phase and quadrature components [6]; we have also looked into the 
Karhunen-Loeve expansion as discussed in the previous reports. The 
Walsh transform of EEG signals offers potentials not only in its efficiency 
of representation but also in its enhancement of certain desirable 
features of EEG signals. The Walsh transform is a representation in 
terms of a set of Walsh functions which were originally developed by 
J. Walsh in 1923 [s] and gained populatiry since 1968 [?] . Preliminary 
results showed that the Walsh spectra of EEG in the various sequency 
bands gave a better indication of the alertness states. 

The representation of a signal x(t) in terms of the set of Walsh 
functions is given by 

\ , 0 ^ t ^ T 

k=0 

where k is the number of zero crossings of Walsh functions. Lackey 
[8] presented a good review of the construction of Walsh functions from 
products of the incomplete, but simpler set of Radamacher functions. 

The Walsh coefficients are obtained by 
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(k, • 

For a digital (or discrete-time) signal which consists of. N 
sampled data points, we may represent the elements of the digital signal 
by 


'o 


N-1 


x(m) = E wal (m,k) , m = 0, 1, 2, N-1 


The Walsh coefficients are then given by 


N-1 

x(m) wal (ni,k) , k = 0, 1, 2, N-1 

m“0 

The N-tuple formed by Walsh coefficients in ascending order of the 
index k is called the sequency ordered discrete Walsh transfonn; 
i.e., X = •*-> ‘ Using the property that the discrete 

Walsh function is symmetric with respect to its arguments; i.e., 
wal (k,m) = wal (m,k) , we can write the inverse Walsh transform as 

N-1 

x(k) = ^ ^ f(m) wal (k,m) , k= 0, 1, 2, ... , N-1 

m=0 
t h 

where x(k) is the k element of x . It is thus shown that the 
forward transform and the inverse transform have exactly the same 
operations. The Walsh transforms are implemented through the fast 
Walsh transform (FWT) algorithms. 

Typically, FWT schemes exploit the stricture of the FFT except 
that the role of 

2Ti 

^ N 

W ^ e 
N 

in the FFT [9] is replaced by -1 in the FWT. This makes FWT so 
much faster than the FFIE The FWT algorithm requires only addition 
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and subtraction. Notable algorithm for the computation of FWT has been 
described in [lO] . 

There are some important differences between signals represented 
in the Walsh domain and those in the Fourier domain. First of all, 
the digital signal N-vector x is represented exactly by the first 
N Walsh functions. This means that sampled data of finite duration is 
also sequency limited. This is, of course, not true in the frequency- 
domain representation where a time'-limited signal is not band limited. 

One can obtain finite segments of a signal x(t) by multiplying 
the signal by a ‘'window function^’ w(t); that is, - x(t) • w(t) , 


where 


1, 0 ^ t ^ T 


w(t) == 


Ko th^^t 



Now, let the Fourier transform of x(t) be X(f) ; i.e., 


F[x(t)] == X(f), then F[x(t) • w(t)] = X(f) * W(f) 


where 


W(f) = K 


sin af 
af 


and K and a are constants; and * denotes convolution. This result 
shows that the Fourier transform is distorted through the convolution 
with the frequency-domain window function. It is well known that the 
error may be minimized by applying a corrective weighting function 
such as Hamming window. This extra processing is used solely to correct 
the distortion effects of the data windows. Such distortion does not 
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occur in the Walsh domain because the window function is one of the basis 
functions; namely^ w(t) == wal (0,t) . This is easily seen by the 
following reasoning: Taking the same data window w(t), we obtain for the 
finite segment of the signal x(t) as 

CO 

x(t) • w(t) - ^ ^ wal (k,t) * wal (0,t) 
k=0 

Since the complete set of Walsh functions is closed in multiplication 
and wal(k,t) * wal (m,t) - wal (k(^m,t) where denotes module 2 

addition with no carry; i*e^, 

10 0.1 
101=0 
0 (?) 0 = 0 , 

wal (k,t) ■ wal (0,t) = wal (k,t) and hence x(t) w(t) = wal (k.t). 

k 

It is clear that the window function w(t) introduces no distortion at 
all in the Walsh domain^ 

In summary, the Walsh-domain representation has two distinct 

advantages over the Fourier-domain representation for discrete signals 

in addition to its computational speed. First, the data is exactly 

represented by a finite set of Walsh functions, and secondly, no 

distortion is introduced by the finite window function. Therefore, no 

correction such as the use of Hamiming or Hanning window is necessary. 

We shall now turn our discussion to some properties of Walsh functions 

which we may utilize to our advantage. Let us consider a signal 

which is periodic with period 2T and a signal 1’ (t) which is the 

same as except that it has period T as shown in Figure 9. 

In other words, is a time-compressed version of f (t) ; i.e., 

1 
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f^Ct) = f^(2t) , t ^ 0 • 




Figure 9 Time compression of a periodic signal 


Representing the signals f (t) and f^(t) in the Walsh domain 

Jl ^ 


obtain 


and 


fl(t) 


‘ ay 

^ ^ wal (k,t) 
k=0 


* 2 <‘> 



wal 

k 


(k,t) 


The time compression condition ^ 2 ^^^ ^ 


implies that 


CO CO 

wal (k,t) = w wal (k,2t) . 


From the definition of the Walsh functions, we have wal (k,2t) - 
wal (2k, t) ; hence, 

CO CO I 

^ ^ wal (k,t) - ^ ^ wal (2k, t) . 


k:rO 


k ==0 


Therefore, 


1 

= W 
0 0 

1 

== 0 

wl = w 

2 1 


W 3 = 0 

"4 ■= \ 


w, = w, 


k " V2 


k even , 


- n 


O j^ r 
v«. V 


For further compression in time such as f (t) = f (mt), where 

^ 1 

jC 

m = 2 , = 1, 2, ,,,, we obtain 


W 


1 

k 



is an integer 


0 , otherwise 


In essence, it says that compressing time by a power of 2 leaves the 
Walsh spectrum completely unaffected except for intervening zeros. 

This property may be used effectively for the detection of subharnionics 
in a signal. For example, we show the spectrum of a signal which has 
period T - 1.024 seconds in Figures 10 and 11. Figure 10 shows the 
spectrum without subperiods and Figure 11 shows that with 8 subperiods. 
Notice there are 8 intervening zeros in Figure 11, 


19 




P'lrJ-r'Y 




sr our.MCY 



The convergence property of the Walsh representation compares very 
favorably with the Fourier representation. To demonstrate, we computed 
the percentage loss in the-^signal energy (or the normalized mean squared 
error) of a section of EEG about 1 second in length for both the Fourier 
and Walsh transforms. The results are plotted in Figure 12, It is 
clear that the Walsh transform has a better convergence property than 
the Fourier transform in the lower range. In other words, in the range 
between 0 to 40 terms the Walsh transform requires fewer number of terms 
to achieve the same amount of mean-squared error as that required by 
the Fourier transform. Figure 13 shows a typical Walsh and Fourier 
spectra of a record of EEG signals. Some preliminary results have shown 
that the Walsh transform does enhance certain features of EEG signals. 
Further work is in progress for utilizing these features for EEG 
classification purposes. 
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IV. A Digital Phase-Distortionless Filter Operating in the Frequency 
Domain 


In phase uncertainty model for the EEG signals to extract the 
parameter b as a measurement of alertness, it is critical to filter 
the EEG signal by a phase-distortionless filter. Such a filter operating 
in the time domain was used in previous studies. Since this filtering 
requires that the EEG signal convolves with the finite duration impulse 
response of the filter, it is extremely slow. This, of course, prohibits 
real-time usage. A phase-distortionless filter must have linear-phase 
characteristics in the frequency domain. In order for a ditial filter 
to possess linear phase characteristics, the impulse response of its 
system function must satisfy the condition 


h(n) = h(N - 1 - n) . 

Tnis requirement is easily met in digital filtering. It results in a 
linear phase shift which corresponds to a delay of (N-l)/2 sample 
points. The input data x(n) is first partitioned into data segments 
Xj^(n) of length M ; i.e., 


X, (n) 
k 


x(n) 

.0 


kM s n ^ (k+l)M-i 
otherwise 


Hence, x(n) is simply the sum of the 


x^(n) 


or 


00 

x(n) - z Xj^(n) 

k=-co 

and the filtered output y(n) is the result of convolution of the two 


sequences; i.e.. 


y(n) - x(n) * h(n) 


CO 

^^Xj^(n) * h(n) 

K=r~CO 
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Each of the terms in the convolution sum is of length (N+M-1) 

where 

N = filter length 
M = data length . 

For each output point, there are numerous multiplications and additions. 

The astronomical number of such operations renders this scheme extremely 

slow. In order to speed up the process, other means must be found. 

We decided to use the frequency-domain approach. From the above discussion, 

we found that to effect the convolution in the frequency domain, the 

transfoi’ms must be computed on a basis of at least (N+M-1) points. 

The (N-1) and (M-1) locations subsequent to the N filter coefficient and 
th 

the M data point, are respectively filled with zeros before each 

transform is taken. Since the FFT results in a circular convolution 

of the two sequences, the desired output sequence is obtained by reconstructing 

the filtered sections in such a way that convolution is achieved. The 

filtered sections are simply the results of inverse transform of the 

sequence obtained by multiplying the transformed data sequence, term 

by term, with the transformed filter coefficient sequence. The procedure 

for reconstructing the filtered sections is referred to as the overlap- 

add method. This complete operation involves FIT and multiplications 

in one scope for obtaining the whole section of the filtered results; whereas 

the time-domain operation requires multiplications and additions for 

each ouput. point. 

Moreover, additional savings in time could be achieved by this frequency- 

domain operated filtering procedure if both of the signal and filter 

sequences are real. In general, the filtered sections y (n) are 

k 
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Our filter is real; hence, 



and becomes 

Yj^(n) = Re 


^Xj^(n)J * Re|^h(n)j 


+ jjim Xj^(n) * Re 


|^h(n)J I . 


If the data sequence is also real, two separate convolutions with the 
same kernel sequence can be performed simultaneously. Therefore, when 
convolving real signals with real impulse responses, the FFT is capable 
of double performance. This is achieved by loading the data in both 
the real and imaginary parts of the data sequence to be transformed 
and processed with the same impulse response of the filter. Thus, 
we are able to save additional time in the filtering process. The 
double-performance overlap-add method is depicted in Figure 14 in which 
section a is discarded, as it constitutes the inherent delay. The 
subsequent sections in Figure 14 are dealt with in the following manner; 
Unprimed and primed sections are added together term by term and the 


continuous filtered output is obtained by stringing together the 
resultant sections^ 


A few words about the judicious choice of the data length M and 
the filter impulse response length N are in order. There are two 
major considerations; (a) N must be sufficiently large so that the 
ratio of the Hamming weighted impulse response side-lobe width to the 




Figure 14 Double-performance overlap-add method 

main-lobe width is sufficiently small in order to achieve a shaip 
cutoff in the frequency response characteristic; (b) The overlap-add 
procedure described previously will be more efficient computationally 
when N = M 'f- 1 and N + M - 1 = 2^ , where i is an integer. In 
our application, we chose N = 4097 and M = 4096 . This choice does 
satisfy the above considerations. The frequency-domain operated phase 
distortionless digital filter which we have designed and used has 
a passband of 8 to 13 Hz. Its amplitude response characteristic is 
shown in Figure 15* Notice the sharp cutoff feature of this filter. 

We show the flowchart of this filter in Figure 16, and the resultant 
filtered output of a section of EEG signal by the use of frequency- 
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Amplitude 





Figure 15 The amplitude response 
of the phase-distortionless bandpass 
digital filter for the a-band 
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Figure 16 Flo^vchart of the phase-distortionless 
digital filter 


























domain filtering- in Figure 18. From these figures we observe that the 
resultant outputs are identical. 

From previous discussion, we may conclude that the phase~distortionless 
filtering in the frequency domain is versatile and faster than the time- 
domain filtering. The use of FIT permits the filtering to be accomplished 
on real-time basis; whereas the time-domain filtering for the same 
length of data will take at least fifty times longer to accomplish. 

In other words, the frequency-domain filtering needs only l/so of the 
time needed for the time-domain filtering. Besides, the FFT is also 
capable of double performance discussed previously. We may use this 
double-duty feature to filter simultaneously either two independent 
channels of EEG signals or one channel of EEG signals in two different 
bands of interest. 
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V. 


Conclusion 


The work has been progressing along the lines as proposed. Our 
preliminary results presented in this report show great promise. The 
tests for stationarity in EEG signals are designed for the determination 
of the length of signal interval in which spectral analysis is valid and 
have great potential in extracting the dynamic characteristics of the 
Signal process. This will, undoubtedly, help in the detection of 
Changing states of alertness. In the area of searching for new represen- 
tation, we described the use of Walsh transform and showed its advantages. 
In the modeling of EEG signals, we described, in detail the phase- 
distortionless digital filter in the frequency domain. Its versatility 
and speed have been proven from our results. This filter can be operated 
in real time in contrast to the time-domain filtering techniqued used 
previously which is fifty times slower. Without question, this filtering 

technique will speed up the estimate of the parameter b in the model 
for alertness measure. 


34 



References 


Kawabata, N., nonstationary analysis of the electroencephalogram 
IEEE Trans . , Vol. BME--20, No. 6, pp . 444-452, November 1973. 

Anliker, J. E., Friedlander, S., Mansfield, F. , and Finger, H. , 
^^Nonstationary processes in the electroencephalogram during sleep”, 
Proc. of the 27th Annual Conference on Engineering in Medicine and 
Biology, p. 145, October 1974. 

Walsh, J. L. , ”A closed set of orthogonal functions”, Am. J. > Math. , 
Vol, 55, pp. 5~24, January, 1923. 

Floyd, R. V., Lai, D, C., and Anliker, J. E., ”A model for the 
photically stimulated electroencephalograhic signals", Proc. of 
the San Diego Biomedical Symposium , Vol. 12, pp. 5--16, 1973. 

Lux, R., "A statistical representation of the alpha rhythm in the 
human EEG with applications to the classification and prediction 
of behavior states", Ph.D, Thesis, University of Vermont, May, 1974 

Ein-Gal, M. and Lai, D. C., "Real-time EEG analysis and monitoring 
using in-phase and quadrature components", Proc. of the 26th 
Annual Conference in Medicine and Biology , p. 401, 1973. 

Haimiuth, H. F. , "A generalized concept of frequency and some 
applications", IEEE Transactions on Information Theory , Vol. IT-14, 
No, 3, pp. 375«382, May 1968. 

Lackey, R. B., and Meltzer, D., "A simplified definition of Walsh 
functions", IEEE Transactions on Computers , pp. 211-213, February, 
1971. 

Ulman, L. J., "Computation of the Hadamard transform and the R 
transform in ordered form", IEEE Transactions on Computers , pp. 359- 
360, April, 1970. 

Shanks, J. L. , "Computation of the fast Walsh-Fourier transformV, 
IEEE Transactions on Computers, pp. 457-459, May, 1969. 



